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ارزیابی آثرات توپ و گرافی بر تفرق امواج لرزه‌ای بااستفاده از روش المان مرزی* 
(یادداشت پژوهشی) 


محسن ایثاری(٩‏ رضا تاری‌نژاد!" 


چکیده در این تحقیق برای ارزیابی اثرات توپوگرافی از روش المان مرزی سه‌بعدی که یکی از فنون عددی ب سیار قدرتم ند برای تحد یل 
مسائل به‌ویژه موضوعات مربوط به انتشار امواج است, استفاده می‌شود. نتایج به‌دس تآمده از تحد یل تو پوگرافی: های مد لف .به اشکل: های 
مثلثی, ذوزنقه‌ای, نیم‌دایره‌ای و نیم‌بیضوی تحت اثر امواج لرزه‌ای 511 نشان می‌دهد که یکی از مهم ترین عوا مل تأثیرٌ گذار بر الا گوی رف تار 
دینامیکی دره پارامتر هندسی (شکل و عمق) می‌باشد, که | ین عوا مل در فر کانس های مد لف | ثرات مت فاوتی از خود ذ شان «می‌ده ند. 
علاوه‌براین؛ اثر پارامترهای دیگر مانند زاویهةٌ برنخحورد و همچنین خحواص مصالح دره (مدول لاستیسیته, میرا یی و ضریب پوا سون) ذ یز مورد 
ارزیابی قرا رگرفتند. 


واژه‌های کلیدی توپوگرافی المان مرزی, تفرق امواج. اثرات ساختگاه هندسة دره. 
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*تاریخ دریافت مقاله ٩۳/۱۲/۲۹‏ تاریخ پذیرش آن ۹۵/۲/۱۵ می‌باشد. 
(۱) کارشناسی ارشد دانشکده مهندسی عمران دانشگاه تبریز. 


(۲) نویسنده مسئول: استادیان دانشکده مهندسی عمران دانشگاه تبریز. که مزرجام) ۵ 4هزممتیها ۲ : اتقطرظ 


۹ ارزیابی اثرات توپوگرافی بر تفرق امواج لرزه‌ای بااستفاده از روش المان مرزی 


مقد مه 
مطالعات زیادی درخصوص پديدهة تشدید توپوگرافی 
روی امواج لرزه‌ای توسط محققان مختلف انجام گرفته 
است. یکی از نمونه‌های بسیار واضح درخصوص این 
پدیده. گسیختگی‌های سنگی به‌وجودآمده در تکیه‌گاه 
سد پاکویما واقع در کالیفرنیا در طول زمین‌لرزه‌های 
بزرگ می‌باشد. بررسی‌های صورت‌گرفته نشان می‌دهند 
این گسیختگی‌های سنگی مربوط به پدیدهة اثر 
توپوگرافی می‌باشد. رکوردهای زمین‌لرزهٌ ثبت‌شلده در 
زلزل ۱۹۹۶ نورتریج و ۱۹۷۱ سن‌فرناندو و بالای کوه و 
قعر در نزدیک مخزن و سد پاکویما نشان می‌دهند که 
ماکزیمم شتاب زمین در بالای کوهی که جناح چپ سد 
را تشکیل می‌دهد در دو زلزلة مذکور به‌ترتیب برابر 8 
۸ و ع ۱/۰۵ بوده است. در حالی که در قعر دره 
۳۶ و 2 ۰/۸۹ اندازه‌گیری شده است. بدین ترتیب 
بزرگ‌نمایی توپوگرافی برای زمین‌لرزهٌ سال ۱۹۷۱ زلزلة 
سن‌فرناندو و زمین‌لرزه سال ۱۹۹۶ زلزلة نورتریج 
مشخص شدند [2 ,1]. موسسیان و دراوینسکی در سال 
۰ به بررسی مسئلهٌ بزرگ‌نمایی امواج الاستیک در 
یک درهٌ سه‌بعدی به‌شکل دل‌خواه پرداختند و پاسخ‌های 
حالت دائم و گذرا را به‌دست آوردند 4 ,3]. یوان و 
جیانگ به‌کمک گسسته‌سازی به‌شکل نوار محدود 
فرمولاسیون روش المان مرزی در حل مسائل تفرق از 
یک دره سه‌بعدی به‌شکل دل‌خواه را استخراج کردند. با 
ایین روش حجم محاسبات کاهش می‌یابد و دقفت 
محاسبات درمقایسه با روش‌های مرسوم المان محدود و 
المان مرزی بالا می‌رود [5]. اشراقی و دراوینسکی 
توانستند بااستفاده از بسط توابع موج در مختصات 
کروی» مسئله تفرق امواج هارمونیک در برخورد با دره 
سه‌بعدی را مورد تحلیل و ارزیابی قرار دهند [6]. از 
موارد دیگر تحقیقات صورت‌گرفته در اين زمینه می‌توان 
به تحقیق تاری‌نژاد و همکاران اشاره کرد که به بررسی 
و آنالیز تآثیرات بزرگ‌نمایی توپوگرافی بر ساختگاه 
دره‌ای بااستفاده از روش المان مرزی سه‌بعدی 


پرداخته‌اند. آنها با درنظر گرفتن بعضی از پارامترهای 


نشریة مهندسی عمران فردوسی 


تأثیرگذار بر بزرگ‌نمایی جابه‌جایی حاصل از زلزله مانند 
زاویه و فرکانس امواج. به بررسی جداگانة ايين عوامل 
پرداخته‌اند [۷). کمالیان و همکاران مطالعات پارامتریکی 
را به‌منظور بررسی رفتار لرزه‌ای تبه‌های نیم‌سینوسی 
نیمبیضی و ذوزنقه‌ای‌شکل دربرابر امواج قائم 5۷ و ۳ 
انجام دادند. در این مطالعات که بااستفاده از روش المان 
مرزی انجام گردید تأثیرات نسبت شکل. نسبت طول 
موج برخوردی به ابعاد عارضه و ضریب پواسون به‌طور 
جداگانه در پاسخ لرزه‌ای تبه‌های یادشده مورد بررسی 
قرار گرفتند [8]. کمالیان و همکاران در مقاله‌ای تحت 
عنوان آنالیز رفتار دوبعدی ساختگاه در حوزة زمان 
برروی توپبوگرافی سطوح بااستفاده از روش المان 
مرزی, به بیان فرمول‌بندی پیشرفته در حوزهُ زمان؛ 
توسط روش المان مرزی دوبعدی برای حالت 
سوه ینامک مین بدا ان [9ا پکتی آدیگی از 
تحقیقات انجام‌ش ده توسط کمالیان و همکاران 
تحلیل‌های پارامتری متعددی است که روی دره‌ای 
به‌شکل نیم‌سینوسی تحت اثر امواج فشاری و برشی 
منتشرشده در راستای قائم در حالت دوبعدی بااستفاده 
از روش المان مرزی مستقیم در حوزهْ زمان انجام دادند. 
طول موج و هندسة ساختگاه و در درجه‌ای پایین‌تر 
مشخصات مصالح و نوع موج مشخصه‌های حاکم بر 
الکو و کاس ات فروها فتاسانی تنل ۲۱0 

امیدوار و رحیمیان در مقالات به بررسی تفرق سه‌بعدی 
امواج هارمونیک مسطح برروی سطوح با توپوگرافی 
خاص بااستفاده از روش المان مرزی سه‌بعدی پرداختند. 
نشان داده می‌شود که برای بررسی هرچه دقیق‌تر 
حرکات سطحی زمین در توپوگرافی‌های مختلف باید 
تمام پارامترهای موْثر مانند هندسه محل. خصوصیات 
مکانیکی مصالح زمین‌شناسی موجود (دانسیته. نسبت 
پواسون. مدول برشی) نوع امواج» آزیموت و زاوية 
برخورد و همچنین فرکانس موج درنظر گرفته شود 
[11,12]. از دیگر مطالعات صورت‌گرفته در این زمینه 
می‌توان به پژوهش‌های گتمیری و همکاران اشاره کرد 


که بااستفاده از روش المان مرزی بهبودیافته. به بررسی 


سال سی‌ام» شمارة دو. ۱۳۹۲ 


محسن ایثاری- رضا تاری نژاد 


رفتار الاستودینامیک دوبعدی سطوح شیبدار در حوزهُ 
وان پرداتبان, این ریق براساس کر ناهفتون 
انتگرال‌های زمانی است که منجر به کاهش سریع 
تعرآب‌شای اشانی با کلذکت سا مب شود [14 :13 
سهرابی بیدار و همکاران به بیان فرمول‌بندی روش المان 
مرزی سه‌بعدی در حوزه زمان برای تحلیل رفتار 
نراشگاه. پرداستن. رای ارژیانی ذفتته فرمنول شتدی 
ارائه‌شده رفتار لرزه‌ای توپوگرافی‌های مختلف به‌شکل 
دره و تپه که دربرابر برخورد امواج 3,٩‏ قرار گرفت 
ارزیابی شد. نتایج به‌دست‌آمده از تحلیل سه‌بعدی رفتار 
لرزه‌ای نشان داد که تیه سه‌بعدی متقارن دارای پتانسیل 
بزرگ‌نمایی بیشتری درمقایسه با تپه‌های نا متقارن 
سه‌بعدی و دوبعدی می‌باشد [15]. سهرابی بیدار و 
همکاران بااستفاده از روش المان مرزی سه‌بعدی 
مستقیم در حوزه زمان پاسخ درهُ سه‌بعدی گوسی‌شکل 
را به‌صورت پارامتریک مورد مطالعه قرار دادند. طول 
موج و هندسة ساختگاه و تا حدودی نوع موج مهاجم 
پارامترهای موثر بر بزرگ‌نمایی این نوع تپه‌ها شناسایی 
شدند. دره‌های دوبعدی بزرگ‌نمایی و کوچکنمایی 
کوچک‌تری را نسبت به دره‌های سه‌بعدی با هندسه و 
اتعاه مانه تقان بر دص 16 )رش یی سار که 
براساس روش المان مرزی سه‌بعدی در حوزهٌ فرکانس 
می‌باشد مطالعات پارامتری متنوعی روی دره با 
هندسه‌های مختلف انجام می‌گیرد و اثرات پارامترهای 
مختلف در پدیدة تفرق امواج لرزه‌ای مورد بررسی قرار 


می‌گیرند. 


حرکات سطح آزاد زمین و دره 
گام نخست در بررسی رفتار دینامیکی سازه‌های واقع در 
ساختگاه‌های دره‌ای در اثر زلزله» شناسایی امواج 
ایجادشده و محاسبهٌ جابه‌جایی‌های به‌وجودآمده در 
محیط نیم‌بی‌نهایت زمین مسطح (بدون وجود دره) در 
محل تکیه گاه‌ها می‌باشد. در حالت سه‌بعدی فرضص 
می‌شود که بردار عمود بر جبهه امواج لرزه‌ای حجمی 
با بردار انتشار موج؛ زاویه‌های ,8 با محور قائم (2) و :0 


سال سی‌ام» شماره دو» ۱۳۹۹ 


۹6۵ 


بامحور طولی دره 09 می‌سازد. در شکل (۱) 
مشخصات سیستم و یک مقطع عرضی دل‌خواه از 
ساختگاه به‌همراه زوایای برخورد افقی و فائم امواج 


تشان داده شده است. 


شکل ۱ مشخصات سیستم و یک مقطع عرضی دل‌خواه از دره [۷] 


در دستگاه مختصات ۳-۰-7 جابه جایی‌های 
امواج برخوردی هارمونیک در هر نقط؛ دل‌خواه از 
ناحيةٌ نیم‌بی‌نهایت زمین را می‌توان توسط سه مژلفة 
جابه‌جایی زیر معرفی کرد. 
0 انم (امم 

که درآن ۸ دامن جابه جایی» ه فرکانس موح, 
,0عنول ی /ه) < 1 و ,۵/6 ۷ ببسرای ام‌واج 
فشاری و ,9عنعل0/6)- 1 و ,۵/6,(090)< ۷ برای 
امواج برشی می‌باشند. م6 و :6 به‌ترتیب سرعت امواج ۳ 


و ٩‏ فطل ق بردای ۲ عبارخ ات از؛ 


| ,0.6050 ,0 «زم) 
1 ,0 0,510 ,,0 05 ِ ۳ 
0 .01) 


ترا ۱۳۹ 101 
)۳( 6 - ۷ 101 
6 - ۳ 101 


جاب‌جایی‌های میدان آزاد زمین با مجموع 
جابه جایی‌های برخوردی و جابه‌جایی‌های ناشی از 
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۹1 ارزیابی اثرات توپوگرافی بر تفرق امواج لرزه‌ای بااستفاده از روش المان مرزی 


امواج منعکس شده برابر می‌باشد. در دستگاه مختصات 
۷-7 که امواج در صفحهٌ ۰-7 منتقل می‌شوند 
جابه‌جایی‌های حالت آزاد را می‌توان به‌صورت زير بیان 
نمود. به‌منظور اختصار در رابطه‌های ارائه‌شده درادامه. 
ضریب "لو حذف شده است. 
0 
+ یومع و9 صف)ملی 1 (),ن) 
0 


0 
7(۲ 0-6090 طلوت) و1 
۷ ۷ «ع 1 


0 


(۳) 


معادلة حاکم بر حل مسئله 
معادلات حاکم بر مسائل الاستودینامیک از روابط ناویر 
تبعیت می‌کند که به‌صورت برداری و براساس مژلفه‌های 
سرعت امواج به‌شکل زير قابل ارائه است. 


3 کت ۷ 


22 


۲۷ )۲۰۵( ۰-05۲۷۰۲۷۰ 


در این رابطه 0+2 عه و کم )یه 
به ترتیب سرعت امواج فشاری و برشی در محیط ۱ 
نیروهای حجمی در جرم واحد و ا بردار جابه‌جایی‌ها 
می‌باشند. ضمن این‌که 7 و ۸ ثابت‌های لامه و م جرم 
حجمی جسم موردنظر می‌باشند. درصورت استفاده از 
روش تحلیل در حوزه فرکانس, بردار جابه جایی‌ها برای 
تحریک هارمونیک با فرکانس ۵ به فرم زیر درمی‌آید. 
2 ع(17)0 - )2 

در رابطة فوق. " دامنة بردار جابه‌جایی‌ها در 
نحوزم فرکانس. من باشد. بااسفاده از رابطتة (ه) فعادلة 
برداری (4) به شکل مستقل از زمان به‌صورت رابطة (1) 
فان وی نی 
)0( ان 3۷0۲-۰2۷۷۷ 

معادلاً انتگرال مرزی حاکم بر مسئله می‌تواند از 


تلوری تقابل دینامیکی به فرم زیر به‌دست آید. 
۹2 ۳ 5 آذام | + نوی 
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در این رابطه آ». مولفه‌های تغییرمکان در نقطه 
مرزی 1 ا و ۴ مژلفه‌های تغییرمکان و ترکشن روی 
تمام مرز. "" و 2 جواب‌های اساسی تغییرمکان و 
ترکشن روی مرز در اثر بار واحد متمرکز در نقطه 1 
مررباشتد: 

ضریب مستقل » معروف به ترم پرش وابسته به 
هندسهٌ خاص مرز در نقط؛ 1 می‌باشد و می‌تواند از 
ترکشن صفر جسم صلب محاسبه گردد [18 ,17]. 

زمانی که مرز به تعداد 06 المان گسسته‌سازی شد 
با جای‌گذاری معادلات پارامترهای گسسته‌سازی‌شده در 
معادلهٌ (۷) معادله زیر به‌دست خواهد آمد: 
۸( «ا یل < ِ بطم 5 زا + لول 

نا ور از 

که ,۳ نشان‌دهندة سطح المان ژ می‌باشد. معادلة 

فوق را می‌توان به فرم زیر بازنویسی کرد. 


9 + ۵۳ - ٩ نوئی‎ 9 
۳1 


۳۱-1 
بنابراین برای تمام گره‌های 1می‌توان سیستم 
معادلات را به فرم کلی زیر بیان کرد. 
۱۰ نحل ۲۱ 


جواب‌های اساسی 

جواب‌های معادلة (1) به‌ازای بار نقطه‌ای هارمونیک با 
دامن واحد که درجهت اختیاری (بردار واحد) اعمال 
گردد. جواب‌های اساسی يا توابع گرین نامیده می‌شوند. 
این توابع بااستفاده از تجزيةٌ هلمهولتز معادلات حاکم 
به‌دست می‌آیند. درادامه جواب‌های اساسی مربوط به 
تغیبرمکان و ترکشن ارائه می‌شوند. این جواب‌ها در 
منابع موجود همواره با غلط‌های تایپی همراه بوده است 
که در اینجا فرمول‌های دقیق آنها بااستفاده از محاسباتی 
که نگارنده انجام داده است ارائه می‌گردند. 

(۱۱ ار ب سا < ,1 


که در آن ۷ و بر برابر مقادیر زیر می باشند. 
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محسن ایثاری- رضا تاری نژاد 


و[ " 1 


احض 2 اج 
ره ۳ 3 " 3 ره (تر ۲ 3 3 


7 یج) < 
1 تا تا 1 تیا 1 ِ 


همچنین ۰-4 و نز نشان‌دهندة تابع دلتای 

کرونکر و - با : به‌تریب اعداد امواج 
1 

قارف پرشی ماه آبتش‌هتانی که ابا ,کاب 
می‌باشند به‌منظور مشتق‌گیری نسبت به جهت موردنظر 
می‌باشد. لازم به ذکر است از روی جواب اساسی 
مربوط به تغییرمکان می‌توان جواب اساسی مربوط به 
ترکشن را محاسبه نمود. این کار بااستفاده از روابط 
موجود بین تنش؛ کرنش و تغییرمکان و نیز رابطٌ تنش 
و ترکشن (قانون استوکس) صورت می‌گیرد. جواب 
اساسی مربوط به ترکشن محاسبه‌شده از روی جواب 


اساسی تغییرمکان به‌صورت زیر می‌باشد [17-18]. 


0 2 0 1 
م2 برر 22 - ر, )2 ّ (ر1 ۳ 2 


0 7 . باه 7 ۵ 17 
۲ 2- 2) + بط بت 2 
ارت( ۲ 4۲ ۷ ت ۱۹ 


4 


که در آن 0۷/۵۲ و 4/۵ به‌شکل زیر می‌باشند. 


(لرئت)صه 3 3 2 1 
/ ۱ 1 ۲ 
1 ی طر1 1 ك 
رل 3 3 0 ۰ 
1 1 ۱ 1۳ ۱ بط 
۱ 9 4 0 
۱ ۳2 شب ۱ ,1 ۳ 1 
1 تیک رک 1 و61 
2 
> 9 9 4 
دس سب یب + + ,۲16 
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برآورد انتگرال‌های تکین 
برای انتگرال گیری از هسته با تکین بودن از نوع ضعیف 
مشابه آنچه در هسته تغییرمکان وجود دارد بااستفاده از 
تکنیک تبدیل انتگرال‌گیری روی المان مستطیلی به 
انتگرال‌گیری روی المان مثلثی می‌توان از تکین بودن 
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سس ۳ (لیعلت) وه 1 ِ ِ ِ یات سك 


۹۷ 


رهایی یافت و از روش گوس- لژاندر انتگرال‌های اجزا 
مثلثی را برآورد نمود. به‌این‌منظور ابتدا المان مستطیلی 
پزسب :لین که بقطه متیع متطیق بر کذام یسک از تقاط 
الماق مرزی باشل به تشاد متاسی السان شاکی تقنشیم 
ده یی( تا مات سای گر 
روی المان مستطیلی به انتگرال‌گیری روی المان‌های 
مثلثی تبدیل می‌گردند. اساس روش تعریف سیستم 
مختصاتی است که جاکوبین در نقطةّ منبع صفر گردد 
به‌صورتی که با صفر موجود در مخرج هسته انتگرال 
رفع ابهام شود. برای این منظور سیستم مختصات مطابق 
شکل (۲) طوری تعریف می‌شود که جاکوبین در نقطة 


منبع صفر گردد. 


8 


رو.وت 


10.0 20.0) 


شکل ۲ 2 و) المان مستطیلی ٩‏ گرهی اصلی و مادر. ع) المان 
تقسیم‌شده به زیرالمان‌های مثلثی با تکینگی در گره شمار؛ُ ۸ 4 و 
6 تبدیل‌های لازم برای انتگرال گیری از المان مثلثی [۷] 


رابطة بین مختصات 5 و ره و مختصات رع « رع 
به‌صورت زیر می‌باشند. 
۱۳( تور + گ(ره-1) ,5+ ع(رو-1) < ع 
که در رابطءة فوق 11,2 و اعداد سللاغ, ع 


نشریة مهندسی عمران فردوسی 


۹۸ ارزیابی اثرات توپوگرافی بر تفرق امواج لرزه‌ای بااستفاده از روش المان مرزی 


نشان‌دهندة گره‌های گوشه مثلث می‌باشند. 
(۱ ,2۸6,05 < وفل,5|ر ل| < ملع 


که در آن ۸ مساحت مثلث در مختصات 
بی‌بعد وچ - رچّ می‌باشد. 

به‌منظور استفاده از روش انتگرال گیری گوس 
حوزه مربعی در سیستم رورو بایستی به مربعی به طول 
وجوه معادل ۲ انتقال یابد. بنابراین بایستی انتقال دوم 
به‌شکل زیر انجام گیرد. 


00۱ گت و 


۱ 2 


بن‌ابراین رابطه ال اه < وال ,| یز| < وعلردل 
به‌دست می‌آید. انتگرال مستطیلی بااستفاده از تبدیل‌های 
فوق به انتگرال‌های مثلشی در حوزة با- با تبدیل 
می‌شود و بااستفاده از روش استاندارد گوس محاسبه 
می‌گردد [18 ,17]. 
بالگ (1 + ۱1[ -000 


1-1 - هاعمهزهآ ,مد 


۱۹) 


صحت‌سنجی کد کامپیوتری به‌کار رفته در تحقیق 

به‌منظور انجام محاسبات لازم در روی ساختگاه‌های 

دره‌ای کد کامپیوتری به زبان فرترن (1۳۸50) نوشته 
شده است. این کد قابلیست حل انواع مسائل 
الاستودینامیک سه‌بعدی را که شامل برخورد امواج 

خی زمیت ال رنه سیاشگاه‌های وودای انیت دازاسخ: 

پسظرن اطیان از ممت‌و آرژبانی کی 186 بران 

محاسبٌ جابه‌جایی‌های میدان آزاد سطح دره یک‌سری 
مقایسه با جواب‌های به‌دست‌آمده در مقالات مختلف 

انجام شده است. 

۱. تحقیق انجام‌شده توسط ژانگ و چوپرا [19] در 
زمینة تحلیل سه‌بعدی تغییرات مکانی جابه‌جایی در 
سطح دره با مقطع نیم‌دایره‌ای منشوری تحت اثر 
موج ورودی ٩۷‏ با زاوية افقی ٩۰‏ درجه و زاويه 
قائم صفر درجه و فرکانس بی‌بعد واحد درنظر 
گرفته شده است (۱. در شکل (۳) مقایسه بین نتایج 


نشریة مهندسی عمران فردوسی 


به‌دست‌آمده از تحلیل عددی با نتایج تحلیل‌های 
چوپرا انجام شده است. روش استفاده‌شده توسط 
جوپر» روش المان مرزی دوبعدی است که با 
انتگرال گیری تحلیلی در بعد سوم برای دره‌ای به 
طول بی‌نهایت به حالت سه‌بعدی تعمیم یافته است. 
ان هی فان ابیت کنو انم ی موی اسان 
مرزی سه‌بعدی استفاده شده است. اختلاف اندک 
بین نتایج را (کمتر از ده درصد) می‌توان به اختلاف 
در روش به‌کاررفته و استفاده از المان‌های دوبعدی 


توسط چوپرا و سه‌بعدی در تحقیق حاضر نسبت داد. 


2 8 2۱209 ۱ ۸ 
0۳۵۵۲۵۱ 8 2۱209 -عل۱ا 9 
۱ 
/ -ع] -.--- 


عا آمطط۸ صعصومعاموزط 
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شکل ۳ مقایسه نتایج به‌دست آمده از تحلیل عددی و نتایج 
ارائه‌شده توسط ژانگ و چوپرا [19] برای موج ورودی 5۷ با 
0 2 ,0 20 ,0 و فرکانس بی‌بعد واحد برای درةٌ نیم‌دایرة 


نشوری 


۲. مقايسةٌ نتایج بدتینای امه از تحلیل عددی و نتایج 
ارائه‌شده توسط کمالیان و همکاران [10] برای موج 
ورودی 5۷ با 0,90 و 0,20 به‌ازای فرکانس 
بی‌بعد مساوی نیم در شکل (4) ارائه شده است. 
نتایج به‌دست‌آمده از تحلیل حاضر با نتایج ارائه‌شده 
توسط کمالیان و همکاران اندک اختلافی نشان 
نس تفس ضستن آسن کنووشن کل کاماه بکساق 
می‌باشد. نتایج ارانه‌شده توسط کمالیان برای درة 
دوبعدی بوده. درحالیکه در تحلیل حاضر مدل 
سه‌بعدی از دره به‌طول محدود (۵ برابر شعاع دره) 


با روش المان مرزی به‌کارگرفته شده است. به‌طور 
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محسن ایثاری- رضا تاری نژاد 


کلی جواب‌های به‌دست‌آمده از مدل‌های دوبعدی 
نسبت به مدل‌های سه‌بعدی اندکی کوچک‌تر می‌باشد 
و این مورد درمقايسة نتایج به‌دست‌آمده به‌وضوح 
نمایان است. براساس مقايسة نتایج انجام گرفته با 
رفرنس‌های متعدد کد المان مرزی توسعه‌یافته در این 
تحفیق صحت سنجی تن و از نتایج با تست | متام از 
آن برای انجام تحلیل‌های بعدی اطمینان حاصل 
گردید. 


۷ ۱ [) سس 
2( ۰۰۰۰۰۰۰۰۰ 
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شکل ۶ مقايسة نتایج به‌دست‌آمده از تحلیل عددی و نتایج 
ارائه‌شده توسط کمالیان و همکاران [10] برای موج ورودی 5۷ با 


0 2 ,0,0۰0 و فرکانس بی‌بعد نیم برای درهُ نیم‌دایرة منشوری 


تحلیل‌های عددی 

در این مطالعه بااستفاده از روش المان مرزی سه‌بعدی 
توسعه‌یافته به بررسی پارامترهای اثرگذار بر پاسخ 
دینامیکی ناهمواری‌های مختلف پرداخته می‌شود. از 
عوامل تأثیرگذار در این رفتارها می‌توان به پارامترهایی 
نظیر هندسه توپوگرافی. جنس مصالح. پارامترهای امواج 
برخوردی (زاویه و فرکانس) اشاره کرد. برای بررسی 
پارامتر هندسهٌ تویوگرافی از چهار دره به شکل‌های 
مثلثی. ذوزنقه‌ای. نیم‌دایره و بیضوی استفاده می‌شود. 


ارزیابی اثرات هند سه ومشخصه‌های امواج بر درة 
به منظور ارزیابی اثر شکل و عمق دره روی تفرق امواج 
از دره‌های سه‌بعدی تحلیل‌هایی روی دره‌هایی با 
شکل‌های مثلثی و درادامه روی دره‌های ذوزنقه‌ای؛ 
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۹۹ 


نیم‌دایره‌ای و نیم‌بیضوی انجام شده است. ابتدا دره‌ای 
به‌شکل مقلت با ننبت عرش قاعده (0) به عمق (ه) 
پرابر ۲۰۲و برای نشان دادن اثرات عمق و بازشد کی 
دره که از فاکتورهای مهم در طراحی سدها به‌شمار 
می‌آید» مدل می‌شود. شکل (۵) نشان‌دهنده دره مثلشی 
مدل‌شده می‌باشد. این مدل تحت برخورد موج برشی 
1 با زاوب؛ة انقی و قائم 40 درجه قرار می‌گیرد. 
به منظور ارزیابی اثرات فرکانس موج زمین‌لرزه این 
تحلیل‌ها به‌ازای فرکانس‌های بی‌بعد مختلف انجام شده 
و نتایج آن در شکل (1) ارائه شده اسست. نتایج نشان 
می‌دهد که با کاهش این نسبت اثرات توپوگرافی برروی 
رفتار دینامیکی در مثلشی افزایش می‌یابد. همچنین 
تراث عمق: و شکل ذره تایعی او فرکانس می‌ناشد: سا 
فزایش فرکانس به‌دلیل کاهش طول موج برخوردی» 
ثرات ساختگاه افزایش می‌بابد. در فرکانس‌های پایین 
(1>0.1) به‌دلیل بزرگی قابل توجه طول موج نسبت به 
بعاد عارضة طبیعی تفرق امواج ناشی از وجود عارضه 


به حداقل می‌رسد و درواقع موج ورودی اثر وجودی 
دره را احساس نمی کند و پدیده تفرق روی امواج قابل 
چشم‌پوشی است. به‌عنوان یک قاعده کلی درصورتی‌که 
ابعاد دره با طول موج امواج برخوردی به آن قابل 
مقایسه باشد تفرق و آشفتگی ناشی از وجود دره روی 
الگوی تغییر مکان‌های به‌دست‌آمده در سطح آزاد دره 
تس مس ی 3 

درادامه به‌منظور بررسی تأثیر تغییرات زاوية موج 
برخوردی به دره روی نتایج به‌دست آمده تحلیل‌هایی 
با زاویهٌ برخورد مختلف انجام شده است و نتایج آن در 
شکل (۷) ارائه شده است. نتایج نشان می‌دهد که تغییر 
زاویه تأثیر زیادی روی الگوی رفتار دینامیکی و 
همچنین مقادیر تغییرمکان‌های به‌دست‌آمده در سطح 
آزاد دره می‌گذارد. برخورد امواج درجهت قائم منجر به 
جابه جایی‌های متقارن حول محور تقارن دره و برخورد 
با زاویه‌های غیر از قائم منجر به جابه‌جایی‌های نامتقارن 
حول محور تقارن دره می کون در شرایط برخورد غیر 
قائم غالبا یک سمت دره که در سمت برخورد امواج 


نشریة مهندسی عمران فردوسی 


۱۰۰ ارزیابی اثرات توپوگرافی بر تفرق امواج لرزه‌ای بااستفاده از روش المان مرزی 


جاب‌جایی‌های کوچک‌تری را درمقایسه با 
جابه‌جایی‌های حالت قائم نشان می‌دهد. 
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شکل 1 مقايسة نمودارهای به‌دست‌آمده از تحلیل‌های انجام‌شده 
به‌ازای فرکانس‌های بی‌بعد ۵/۱ ۰/۰ و ۱ و شیب دیواره ۶۵ ۲۹/۱ 


و ۱۶/۶ درجه 


نشریهةٌ مهندسی عمران فردوسی 
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شکل ۷ مقایس نمودارهای به‌دست‌آمده از تحلیل‌های انجام‌شده 


در زاویه‌های مختلف در فرکانس بی‌بعد واحد 


ارزیابی اثرات هندسة درة ذوزنقه‌ای 
در ذوزنقه‌ای مطابق شکل (۸) به‌صورت سه‌بعدی 
منشوری درنظر گرفته شده است. در این تحلیل‌ها سه 
نسبت مختلف قاعده بزرگ (0) به عمق (۴) تحت اثر 
برخورد موج ٩۳‏ با زاویه ۶۵ درجه و دو فرکانس 
بی‌بعد واحد ونیم مورد توحجه بوده است. نتایج 
به‌دست آمده از این تحلیل‌ها در شکل )٩(‏ ارائه شده 
آسنت. نتایج نشان می‌دهند کاهش نسبت (0/0) 
همان‌طوری که در درهُ مثلشی گفته شد. باعث بیشتر 
شدن آشفتگی‌های ناشی از اثرات توپوگرافی می‌شود و 
میزان این افزایش به فرکانس وارده وایسته می‌باشد. با 
افزایش فرکانس بی‌بعد میزان اين افزایش بیشتر می‌گردد 
و با کاهش فرکانس این تغییرات کاهش می‌یابد. 
قاعده بالای ذوزنقه (0) به قاعدهٌ پایین (2) در 


سال سی‌ام» شمارة دو. ۱۳۹۲ 


محسن ایثاری- رضا تاری نژاد 


فرکانس‌های بی‌بعد مختلف تحت زاویه قائم بر رفتار 
دیانیکی و میتی اسان حایی هااینه اه نقسن: با کناهتن 
نسبت دو قاعده. شکل دره به‌سمت مستطیل و با 
افزآیتن. آن به‌سبت: لت یل عی کته این تخیر شک 
منجر به تغییرات بسیار زیاد شیب دیوارُ دره می‌گردد و 
این تغییرات باعث می‌شود تا اثرات شکل دره که به 
فرکانس وارده نیز وابسته است دچار تغییرات زیادی 
شود. درواقع می‌توان گفت که با افزايش شیب دیواره 
میزان آشفتگی‌های ناشی از اثرات ساختگاه افزایش و 


با کاهین این شیب این میران کاهش ی بان نتایج ان 


کون فک ( اراک یه اس انب دس گنهن 
فرکانس‌های بی‌بعد کمتر از نیم با افزایش نسبت دو 
قاعده مقادیر جابه‌جایی‌ها به مقدار بسیار ناچیزی 
افزایش می‌یابد. باتوجه به مقادیر ارائه‌شده در جدول 
(۱) می‌توان گفت که در فرکانس‌های بی‌بعد نیم و 
بزرگ‌تر از نیم و در نسبت‌های قاعدة بزرگ‌تر از ۲/۵ 
تاد ال ابص هس تون ۳ 
به ۱ کاهش می‌یابد بیشترین درصد افزایش جابه‌جایی 


است که وقتی نسبت دو قاعده ذوزنقه از ۲۰ 


در فرکانس نیم ایجاد می‌شود. 


یسب 


۰۱5۵0 ۰۱۵0 40 0 


20 ۹۵0۵ ٩2۵ 200 220 30 3۵ ۸0 ۵ 


450 ۰4۵0 ۰36۵ ۰300 ۰22۵ 0 


شکل ۸ دره ذوزنقه‌ای 


جدول ۱ مقادیر جابه‌جایی‌های به‌دست‌آمده در فرکانس‌های مختلف در نقطهٌ صفر و تحت زاويه قائم 


25 5 10 20 
2.593 | 2.258 | 1.665 3 
2.712 | 2.351 | 1.9499 | 1 
1.637 | 1.266 | 1.203 ۱ 2 
2.589 | 2.253 | 2.033 | 6 
1.775 | 1.321 | 1.137 | 1 
0.914 | 0.998 | 1.054 6 
1.711 | 1.749 | 1.762 | 3 


8/4 1 125 
][-2 1.985 | 2 
]1-1.5 | 1.376 | 7 
] 1-1 2.005 | 5 
2-075 | 2.378 | 1 
]1-0.5 | 3.576 | 4 
]1<-0.25 | 0.841 | 51 
20.1 | 1.642 ۱ 6 


جدول ۲ درصد تغییرات جابه‌جایی در فرکانس‌های مختلف 


015 05 025 01 


20.96 | 237.04 | 22.56 | 89 


سال سی‌ام» شماره دو» ۱۳۹۹ 


60.14 


1:3 


4 | 38.81 | درصد تغییرات 


نشریه مهندسی عمران فردوسی 
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ارزیابی اثرات توپوگرافی بر تفرق امواج لرزه‌ای بااستفاده از روش المان مرزی 
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۱-05 


05 00 05 
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(۷۲/۳) 66هاوننا 


۸۲۲6 269۳6۲۵اهوزه 8مباامم۸ امعطعموامونط ۸۲۲۵ 66۳060۲هاموز۵ 


عم اموز0 


شکل ٩‏ مقايسة نمودارهای به‌دست‌آمده از تحلیل‌های انجام‌شده در فرکانس بی‌بعد نیم و واحد و زاوية 0 درحه 


نشریةٌ مهندسی عمران فردوسی 


سال سی‌ام» شمارة دو. ۱۳۹۲ 


محسن ایثاری- رضا تاری نژاد 
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10/21 سس 
15 سم مس 


عنام 1و ععواموزط]ا 


عینام۸۳ 66860۴ امون9 


05 1 15 2 


0 
0512006 )۷/6( 


ولینا۸۲۲۵ 269۳۳06۳0۸اموزظ۱ 


2 15 1 05 


شکل ۱۰ مقایس نمودارهای به‌دست‌آمده از تحلیل‌های انجام‌شده در نسبت‌های دو قاعدهْ مختلف 


ارزیابی اثرات هندسة درة نیم‌دایره‌ای 
برای ارزیابی اثرات تغییر عمق در در نیم‌دایره‌ای از دو 
دره به شکل‌های نیم‌بیضوی و نیم‌دایره مانند شکل 
(۱۱۷) استفاده می‌شود. در نیم‌بیضوی دارای دهانه‌ای 
برابر دهانهُ دایره ولی عمقی برابر نصف شعاع دره دارد. 
تحلیل‌ها در سه فرکانس مختلف با مشخصات مصالح 


نسبت پواسون ۰/۳ و سرعت برشی ۵ انجام شده 


سال سی‌ام» شمارة دو, ۱۳۹۲ 


است. نتایج به‌دست‌آمده از این نمودارها نشان می‌دهد 
که در یک فرکانس مشخص تغییر شکل و همچنین 
تغیبر عمق تأثیرات محسوس و قابل توجهی را بر پاسخ 
دینامیکی دره می‌گذارد. نتایج نشان می‌دهد که دره 
نیم‌دایره دچار بزرگ‌نمایی و تغییرات بیشتری درمقایسه 
با در نیم‌بیضوی می‌شود. چون عمق در نیم‌بیضوی 
نسبت به در نیم‌دایره‌ای درمقایسه با طول موج برشی 


نشرية مهندسی عمران فردوسی 


۱ 


کوجک‌تسر است لطا در ایسن حالست اثر دره روی 


جابه‌جایی‌های به‌دست‌آمده کمتر خواهد بود. به‌طوریکه 


ارزیابی اثرات توپوگرافی بر تفرق امواج لرزه‌ای بااستفاده از روش المان مرزی 
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خواهد شد. نتایج مربوط به مقايسة دو شکل نیم‌بیضوی 
و نیم‌دایره‌ای در شکل (۱۲ آورده شده است. 


شکل ۱۲ مقایسة نتایج به‌دست‌آمده برای دره‌های نیم‌بیضوی و نیم‌دایره‌ای تحت امواج 8۳1 با زاویهُ برخورد (0 07 ,90 < 0) 


نشریه مهندسی عمران فردوسی 


سال سی‌ام. شمارة دو, ۱۳۹۲ 


محسن ایثاری- رضا تاری نژاد ۱۵ 


مفایسه بین دو درة مثلثی و ذوزنقه‌ای (10,و,0/2-2.5) آورده شده اسنت. درصدهای ارائه‌شده 
مقایسه‌ای بین دو شکل ذوزنقه‌ای و مثلشی با عمق و در این جدول نشان‌دهندة اختلاف ایجادشده براثر تغییر 
دهانهة برابر انجام شده که نمودارهای موردنظر در شکل نسبت 0/2 در در ذوزنقه‌ای‌شکل درمقایسه با در مثلثی 


(۱۳) ارائه شده است: درم ورد نتایج به دست آمده 1 ست که هم دارای عمق و هم دارای دهانةٌ برابر با در 
مقايسه نمودارهای مثلثی و ذوزنقه ای می‌توان گفت که ذوزنقه‌ای می‌باشد. همان‌طور که از جدول مشخص 
دلیل اختلاف این نمودارها به‌خاطر وجود قاعده در ست بیشترین درصد افزایش در نسبت 0/222.5 و 
ذوزنقه می‌باشد. که باعث می‌شود اثر موج روی این فرکانس بی‌بعد نیم می‌باشد, که به‌اندازه؛ُ 7۱/۶۸ درصد 
شکل بیشتر و تغییرمکان آن زیاد شود. علاوه‌براین در افزایش می‌یابد. همچنین می‌توان نتیجه گرفت که برای 


جدول شمارة (۳) نتایج به‌دستآمده از تحلیل‌های نسبت‌های 0/۵<10 رفتار در ذوزنقه‌ای با دره مثلشی 
انجام‌شده مربوط به یک نقطه دل‌خواه رفن استن بررسی (هم‌عمی و دهانه) یکسان می‌باشد و می‌توان از اثرات 
نقطء صفر) در فرکانس‌ها و نسبت دو قاعده تغییرشکل صرف‌نظر کرد. 


جدول ۳ مقایسة نتایج به‌دست‌آمده از تحلیل‌های مختلف روی دره‌های ذوزنقه‌ای و مثلئی در فرکانس‌های مختلف و در نسبت‌های 
9/222.,,0 


۸ 2.0 5 10 
][-1.5 0.065 05 12 
]1-1 228 137 903 
1205 68 215 111 
] 1203 1697 992 3281 
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شکل ۱۳ مقايسةً نتایج به‌دست‌آمده برای دره‌هایی با شکل‌های ذوزنقه‌ای و منلثی تحت امواج 1 با زاویة برخورد 45 ,0 ,45< ,0 و 


رای تب واه 


سال سی‌ام شماره؛ دی ۱۳۹۱ نشریة مهندسی عمران فردوسی 


۱۹ ارزیابی اثرات توپوگرافی بر تفرق امواج لرزه ای با استفاده از روش المان مرزی 


ارزیابی اثرات حز مصالح بر تفرق امواج اژ ات 


دره‌های سه‌بعدی ۳ ده 
بمتظور ارژبایی: فراک متجساگ: مصالج شامل خسر زب 
میرایی. ضریب پواسون و همچنین مدول الاستیسیته 
تحلیل‌های مختلفی با مقادیر متفاوت از این پارامترها 
انجام شد. نتایج به‌دست‌آمده در دو زاوي؛ برخورد 
مشعلف فر شکل‌های (6) و (۱۵ اراشه شده‌اند. با 


عیامط۸ 01ع82ععهاموزه 
َ 


ه‌ 


(۷/۴) 12066وننا 


مقايس نتایج به‌دست‌آمده از نمودارهای مربوط به این 


پارامترها می‌توان فهمید که تعییر ات پارامترهای میرایی 52 2 عل16۵8 - بل 
ظ - عللا1/08 - لاس سس 
7[ 


25 


و پواسون تأثیری بر الگوی تغییرات تغییرمکان‌های 
رن ندارد ولی در مقادیر به‌دست‌آمده تأثیر راز 


1-1 


۲ علاوه‌براین به‌عنوان یک قاعده کلی می‌توان گفت 


که افزايش میرایی منجر به تغییرمکان‌های بزرگ‌تری 


عبامط۸۳ احعحصععهام‌واط 


۱ 
۱ 
۹ 
۱ 
۱ 
۱ 
۱ 
۱ 
۱ 
۱ 


نسبت به حالت بدون میرایی در راستای وجهی می‌شود 


که در معرض برخورد سای قرار می‌گيرد. در مورد تغییر ۳ 10 5 00 0 0 90 
(۷/۳) ع206وزنا 

پارامتر ضریب پواسون نیز می‌توان گفت که با تغییر این 

پارامتر تغییرمکان‌های به‌دست‌آمده دچار تغییر می‌شود. شکل ۱۶ مقايسة نتایج به‌دست‌آمده تحت امواج 5۳1 با زاوية 


با توجه به نتایج حاصل‌شده لازم است که برای به‌دست برخورد (0 ,0 ,245 ,0) و فرکانس بی‌بعد واحد 


آوردن مقادیر دفیق تغییرمکان دره» مقدار دقیق این 
0 > معط - بل 


پارامترها درنظر گرفته شود. درمورد مدول الا سگیستنثه پب صت 0 دس سس 


1 ع 00770 - لا سس 


کشت از قق بارامعی کیک مر الکوعن تغییرمکان‌های 


تفت الم تأثیر دارد و مانند میرایی و ضریب پواسون 


عمط معطصععهامواط 


تعیین مقدار دقیق آن لازم می‌باشد. 
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20 15 10 05 00 05 10- 15- 20 
(۷/۳) ع26اوزنا 


05 00 
(۲/۳) ع9200زنا 


نشرية مهندسی عمران فردوسی سال سی‌ام شمارهُ دی ۱۳۹۲ 


محسن ایثاری- رضا تاری نژاد 


052 2 علب۸/۵8 - پل 
108 - هس مس لس 30 
28 1/08 با اس 


عبنام‌۸ امعططعهمهاموزط 


20 15 10 05 00 5- 10- 15- 20 
(۷/۴) ع0هاوزنا 


شکل ۱۵ مقایسه نتایج ب‌دست‌آمده تحت امواج 71 با زاويةٌ 


برخورد (45 ,0 ,45 < ,0) و فرکانس بی‌بعد واحد 


نتیجه گیری 
در دو دههٌ اخیر روش‌های مرزی محبوبیست فزاینده‌ای 
پیدا کرده‌اند. با چنین الگوریتم‌هایی تنها شرایط مرزها و 
کسسته‌سازی نها برای حل معادلةٌ موج به‌کار می‌رود. 
این روش از معرفی مرزهای موهومی که در روش‌های 
ناحیه‌ای معمولا مورد نیاز است اجتناب می‌کند. در این 
امواج لسرزه‌ای بااستاده از روش المان 


مراخخ 


مرزی سه‌بعدی مورد بررسی قرار گرفت. نتایج 
به‌دستآمده نشان داد که پارامترهای مختلفی مانند 
پارامترهای موج (فرکانس و جهت برخورد) و همچنین 
پارامترهای هندسة توپوگرافی از عوامل تأثیرگذار بر 
الگوی جابه‌جایی دره می‌باشد. همچنین نشان داده شد 
که میزان اثرگذاری شکل دره و عمق دره به فرکانس 
وارده بسیار وابسته می‌باشد. به‌این‌صورت که با افزایش 
فرکانس اثرات افزايش عمق بسیار محسوس خواهد 
بود. مقدار بزرگ‌نمایی و افزايش عمق نیز با یکدیگر 
ارتباط مستقیم دارند. همچنین از نمودارها مشخص 
است که اثرات مصالح (میرایی و پواسون) برروی رفتار 
دینامیکی و الگوی جابه‌جایی غیرقابل‌انکار است و برای 
افزايش دقت در محاسبات لازم است که درنظر گرفته 
شوند. درمورد دره‌های ذوزنقه‌ای نیز می‌توان گفت که 
به‌دلیل وجود قاعده پایین میزان جابه‌جایی بیشتری 
نسبت به در مثلثی هم‌عمق و دهانهٌ برابر با آن ایجاد 
می‌شود. که این موضوع به‌دلیل افزايش سطح 
توپوگرافی موثر در تفرق امواج لرزه‌ای می‌باشد. 


دود , ممتام۱۷ ممامین یمکتصح ما رالمتاهمی اه ون مصتمعوظ ۵ وتورولحمم تهممتامما و۷۷۰ ,وعبلخ ...1 


(2005) ,1603801081 0۶ مابا)تافصا متصصمتلمت رحمتاها تم ووز<[ 


م۵۵00 تاو رممتاما۱۷ مومت مئتص موم ما فصصونا مق ۵۶ ععمممومک معلقتامطاتقط ر.ظقرمنژمک ...2 


1968(۰) ,۱۲۱۲۷۵۲۹۱۷ 0211101012 مطنام 


۰ 1۱1۳260510021 معط[ ج وه ۷۷2۷۵۵ متافهاظ 0۶ ممتاه‌تگتامصم ری رتافصنه:۲ م2 نک روزووموو۱۷۵ ...3 


1990(۰) ,19 ۷۵۱۰ ,]۳۱6 م2 عحظ طاتقظ ر ععممم‌عمک اصمتوصه 1 :22۲2 


لقصمتعصعصزنا معط ه ها ۷۷/2۷۵۵ عتافقاط که ممتاه‌گتامصی و رتاعصا 1۳2 مه بک بتصم[ رصهلوععووم]۷ . .4 


1990(۰) ,۷۵۱۰19 ,501017۲ 200 فص اتف رد عومممومع ملقاو- معا :ایهم .رمال2 ۱۷ 


۲ ۱۷۵۱00 )معا رتملصمظ لهعتارلممش نمی م۸ وک ۱۷۲۰ رعصتامص 2020 .ل .2 رقطمل2 ری ) مطقلا ۷ .5 


1990(۰) ,19 .7۴,۷۵1( 50۱6 290 عصظ طاتوظ همهم ]۲۱2 2 صا ۷۷/2۷۵ 0۲ عصتماهه 


و0 ها ۷۷2۷۵5 تما چم صرح ظ ٩۷,‏ بل متصمصنهز۲! مصعاظ ۵ عصتمااههی و رتکافط ناتک ۲ رتطعمتطو ...6 


0ص عمط اتف ر طمهممموخ ممتعصهمدظ ممتامم۲ ۷۷۵2۷۵ ۸ :5ه۵ رصق [102وصمصرزنال معتط آ. ملتامصصط رفن2 


٩106.17, ۷۵1۰ 18, )1989( 


۷ _ تاری‌نژاده رضا «بارگذاری زلزله برای سازه‌های واقع در ساختگاه‌های دره‌ای»» رسالهٌ دکتری با راهنمایی دکتر محمدتقی احمدی, دانشکدهة 


سال سی‌ام» شماره دو» ۱۳۹۹ 


نشریه مهندسی عمران فردوسی 


۱۸ ارزیابی اثرات توپوگرافی بر تفرق امواج لرزه‌ای باستفاده از روش المان مرزی 


فنی مهندسی. دانشگاه تربیت مدرس, (۱۳۸۲). 
ممتاح‌نآموصد جوم فامع]۲۳ ممفطای یم رطفقطلفم مه رتملاظتصامطمو ویکل۱۷ بتتعلهل با مطعتلمصه۴ 
2003(۰) ,(۳51 19) 58-70 .۳ ,0.2 ر6 ,۷۵۱ مرول و ع۲11 امممتعصمهزنا م19 ۵۶ آعتاصه۳ 
وتوراقصد ععممموم مان اقممتعممصصتل-۲۳۷۵ متقموملمحطا ]مر ریم رتطحتطمو پیک,]۷ رتتقلهل رب مطعتلمصه۴۵ 
,39-45 .۳ ,(5)2 بقاوع آمعنوه7۱0۵۵و0ک .۷ ,۳۱۳۲ظ ها دعتنتامنتوگ منطمدنه 108۵‏ 0۲ 
2032۰ 
ما۷2 2060و-عص1و- جع 21۱ ۵۶ ماو ممتله‌گتامصض بخ رتطاحتطاهی ر.ظ توت بل مصعتامصصفکز 
-871 .۲۳ ,23 ,عطعص ,۱۷۵ .تمحصیااط یمهم ۷۷2۸۷۵۵ )م06تمصا مصتاممهم۵:ظ باله‌تانه ۷ ما 0۲60عزهانا 
2007(۰) ,6887 
0ص «00ظ منصمصیه۱ظ مصقاط ۵۶ ممتاه‌نام‌صه وه فامع]۱۳ بطاممتومع10 ناد.۱۷ رصوتصتطفگ ر.ظ هنن 
2010(۰) ,ول ,۷۷2۵۷۵۵ ۱01266 
٩۷, 220 ۳ ۷۵۹‏ ,9۳۲ منصمصصتها۲ مصهاظ ۵۶ ممتمه ۲۵6-1۲18۵۵۵9100۵ ۰۱ ۱۷۲۰ رصهتحصصتطفکا,.ظ رهبل‌تن 
۸ 0۶0۲۸۷۶ ۲۷۵۸۲/۲۷۵ 0۷ ۵۳ 70۵۳2۷۸ ۸۵۲۸۸ ۷۵۱۱۵۷۶۰ ۵۱۱۵۷12 ۱۷۲۵۱۱۵۷۵۲۵0 ظ1 
2010(۰) ,605-626 ۲۳۰ ,5 .۱0 ,11 ۷۵۱۰ (۲70(۱92۷6] 
صمناه تا تاصهدا۵ .11 ۱۷۲۵۵۵۵0 ۳/۳۲۳۲ظ 21 260نصتانن0 مه اه مامع]۲ مار متفه .۲ رطموتم 20 .ظ بتتتطادت 
ملدتمطتقظ ‏ مه . ومتمصمه نا لتمو ...رد ورملله ۲۷ رحتصمصتلعد. . لقممتومهم‌صصزن-۱۳/۵0 طا. دامع مال 01 
2008(۰) ,646-661 ۲۳۰ ,28 ,۲۱۱2106601108 
٩۷ ۵‏ 0106۲ظ1 0) 0۲۵0عزهاناو وعم910 ۵۶ عقوم متصصوتهگ ‏ ریک فاطعنا ,۷ باعل رظ تاه 
2007(۰) ,060۴060 .)۷۵ ,افص ۸۷۵۸۱۵۲۰ ۰۱۰ ۲۱۶۰ م۸۲02 صمحصعاظ وتحل‌صاوظ ۷۵0 ۵0مرحصا مق ره 
عا اهومزعطم‌صصتل-عمعط 1 موز ۳۸ظ صتمصرمل مصطز 1 رکل. متلوله! م۷ مصهتامصهک یه ,1027ظ وتو 
,1467-1492 ۰( ,79 ,۱۱8۲8 ۱۷۵۱ ۷۱۵/۱۸۵۲ .۰ ۲۲۱ ماگ منطمهع۰108۵ ۵۶ ولورلممخ ۴۵5۵0886 
2009(۰) 
-99120ا02) 3-۲ 01 مومممومک متصصوتهی. رتتوگهل م1۵22 0مصصصصهطم]۷ 200 ,۷۲ رصهئلمصصهک ره رتهل1ظ-201ط0 
,1429-1442 .3۳ ,183 ۲۲ .ل زا0۵0۵ ,۷۷2۵۷۵۵ اصملتمص مصتاممه(۵:ظ بلله:۵ ۷ 0 ۷21165 5۳02۵60 
2010(۰) 
1993(۰) ,۲۱96۷16۲ ,1۱۳۱۷۱۱۵۴۵۵68 مد فاممصصماظ هل‌صتاه‌ظ و.ژ رمعناعمتمو0<[ 
بل صفص ۳۱ متا و تمه لهمافهاظ صا ولمطاه]۱۷ عممصصماط بته)صتمظ وبظنا رومافعظ ر.ملیت مقتامصه] 
1988(۰) ,004-620019-3 :15۳۷ 
کمم۱۷۱۵ میامن مصارنه ۷ الفتادمگ که ولوراممم امطملفمم‌جصزرز-ععط۱. ریک یه بهععمظ) 200 ما ,220 
امن مه مصلهمصتمص معلمتمطتفظ رممهمو]۳۱۵1 فیامعممعممهم۴ ه صا ممرصهت صیمتتصنا ۵ نامه 


1991(۰) ,911-926 .۳ظ ,20 :وعتصهطو(1[ 
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نشریه مهندسی عمران فردوسی سال سی‌ام شماره دو ۱۳۹۲ 


